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Abstract 



We construct an explicit model of gravity trapping domain wall potential where for the first time we 
can study explicitly the graviton wave function fluctuations for any thickness of the domain wall. A 
I concrete form of the potential depends on one parameter < x < ^ , which effectively parameterizes 

OO ■ the thickness of the domain wall with specific limits x — > and x ^ ^ corresponding to the thin 

and the thick wall, respectively. The analysis of the continuum Kaluza Klein fluctuations yields 
, explicit expressions for both the small and large Kaluza Klein energy. We also derive specific 

OO ! explicit conditions in the regime x > 1, for which the fluctuation modes exhibit a resonance 

' behaviour, and which could sizably affect the modifications of the four-dimensional Newton's law 

at distances, which typically are by four orders of magnitude larger than those relevant for Newton's 
^ ■ law modifications of thin walls. 
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I. INTRODUCTION 



The Bogomornyi-Prasad-Sommerfied (BPS) saturated supergravity domain walls [l|] 
which interpolate between two (BPS) anti-deSitter vacua possess the feature that they can 
trap gravity [2]. The original infinitely thin, Zi symmetric domain wall solution 3] in 



five-dimensions can be treated explicitly to obtain the qualitative form of corrections to 
four- dimensional Newton's law [2i]- Subsequent studies extended some further analysis to 
features of the finite thickness walls (see [8), i9|] and references therein) as well as numerous 
other subsequent studies of special models, primarily employing numerical analyses of fixed 
thickness solutions (for recent studies see e.g., [lo| and references therein) (l^ . 

The purpose of this paper is to advance the analysis of the gravity trapping for domain 
walls in an explicit analytic model, where the graviton wave function fiuctuations can be 
studied explicitly for any thickness of the domain wall. The Schrodinger potential depends 
on one parameter < x < ^, which effectively parameterizes the thickness of the domain 
wall. The two limits x — and x ^ | reproduce the infinitely thin and thick wall lim- 
its, respectively. This is the first analytic explicit example where the graviton wave function 
can be obtained explicitly for any value of allowed x, both for small and large values of 
Kaluza Klein masses of fiuctuations. Intriguingly, we also find the resonance behaviour in 
the regime x > 1 and provide an explicit analytic study of resonances which in turn can 
sizably modify the four- dimensional Newton's law at distances that are typically four orders 
of magnitude larger than those relevant for modifications in the thin wall set-up. 

The paper is organized in the following way. In section [III we present the concrete form 
of the gravity potential in the background of the BPS domain with a finite thickness, inter- 
polating (in Z2 invariant way) between two anti-deSitter vacua with negative cosmological 
constant A. We discuss in detail the matching conditions on the metric and the explicit form 
of the potential that depends only on one free parameter < x < |, effectively parameteriz- 
ing the thickness of the wall. In section UTTl we study the graviton wave function fluctuations 
and in particular focus on the explicit form of the probability density at the center of the 
wall. In subsection IIII Al we obtain the analytic expressions for the probability density at 
the center of the wall both for small and large values of Kaluza Klein masses. In subsection 
IIII Bl we analyse quantitatively the conditions under which the wave function exhibits the 
resonances. In section IIVI we study the implications for the modiflcation of Newton's law, 
and in particular obtain the analytic form of Newton's law modiflcations both in thin (x < 1) 
and thick (x > 1) regimes. Conclusions and proposal for further studies are presented in 
section |Vl 



II. FINITE THICKNESS DOMAIN- WALL MODEL 

For the sake of simplicity the model that we shall discuss is chosen to arise from a Z2 
symmetric flnite thickness domain wall that interpolates between two BPS vacua with the 
negative cosmological constant A in flve- dimensions p^]. Our motivation is to present an 
explicit, integrable model of the Schrodinger potential for the graviton wave function modes, 
which would be explicitly parameterized by the thickness of the wall. Such a model would 
in turn allow for an explicit analysis of the graviton fluctuation modes and the subsequent 
analysis of a modiflcation of four- dimensional Newton's law, thus providing a unifying, 
analytic approach which addresses these effects as a function of the wall thickness. 
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We choose the following form of the potential: 

; 
15 



V{z) = -Vo; |z|<f 



V{z) = , \z\>l. (1) 



2 ' 
d 

This potential enters Schrodinger equation for thegraviton wave function (for the derivation 
of the wave equation and other details see, e.g., [91]): 



+ V{z) ipm{z) = m ipmiz) , (2) 



dz^ 

where the graviton wave function hf^iy = rjfj^yip^z) exp(iA;^x^) and k^k^ = m? parameterizes 
the Kaluza Klein energy associated with the four-dimensional momentum A;^. 

The potential V{z) is related to the conformal factor A{z) of the domain-wall metric: 

ds^ = exp[-A(z)] (-de + dz^ + ^dxl\ (3) 



through the following relation (see, e.g., j^): 

Along with the boundary conditions ^(0) = 0, ^(0)' = and A{z) — > 21og(fc|2;|) as 
|z| ^ oo the relation (jl]) determines for the potential ([1]) the following form of the metric: 

A{z) = -^log(cosv^lzl); \z\<l, 

Aiz) = \og[em+P?]: \z\>l, (5) 



where k = y — f (see, e.g., [llj) and the five-dimensional Planck constant Mpianck^ was set 
to 1. 

The continuity of the metric and its derivative at the junction \z\ = | imposes the 
following two conditions among parameters Vq, d and f3: 

k{^+P) = cos(^)-2/^ 

i^+pr = IvvotM^)- (6) 

These two conditions can be viewed as fixing ^ + P and Vq in terms of one free parameter. 

%d 



(7) 

The parameter x has a range {0, |}. The infinitely thin wall corresponds to the following 
limit of parameters: 

3 

x^O; d^O, Vo^ oo, Vod = - = 3k, (8) 



while the infinitely thick wall corresponds to: 

x^f; d^oo, Vo^O, VVoP^-x, {^^P + x)' {^^Y = ^ . (9) 



3 



III. WAVE FUNCTION 



The solution of the Schrodinger equation ([2]) has eigenvalues > 0, with m = 
corresponding to the graviton bound state (see e.g., [9|): 

M^) = No e(-i^(^)) , (10) 
where the normalisation constant Nq is fixed by the following relationship: 
k 

—J = I sin(x) cos(x)~^/^ (x + sin(x) cos(x)) + cos(x)'^/^ . (11) 

In the thin wall limit (x — > 0) and thick wall limit (x ^ |) the normalisation coefficient Nq 
takes the following respective forms: 

k k IT n , . 

— ^^1, as X — : — ^ r^TTT , as x^— . (12) 

' ' 3cos(x)5/3' 2 ^ ^ 

The solution iprn{z) of ([2]) in the continuum, > 0, has the following form: 



ipm{z) = A^cos{K\z\); \z\ < f 

ijrniz) = NmV\z\+P [amY2{m\z\ + h„j2{m\z\ +/?)]; \z\>{. (13) 



Here K = ^/w? + Vq and the coefficients and hm satisfy + 6^ = 1. The normalisation 
constant is determined by employing a regulator \zc\ oo along with the asymptotic 
form of iprn{z) as \z\ oo: 

ijm{z) = Nm^J^[am siYi{m{\z\ + /3) - ^) + 6^ cos(m(|z| +13)- f )] . (14) 

In the limit \zc\ — >■ oo the asymptotic wave function f[T^ ensures the dominant contribution 
to the wave function probability, and thus Nm = \f^'i which in the continuum limit 
becomes = ^/y? i-e- is replaced by vr. 

The matching of the wave function and its first derivative at \z\ = ^ determines the 
coefficients am, &m (ctm + = 1) and A^- In particular, the key expression in determining 
the deviations of Newton's law is the probability density of the wave function at the center 
of the wall: 

l^^(O) |2 = Al = ^ ^ ^ , (15) 

where: 

C = cos(yo), 5* = -ir sin(yo) , (16) 

and 

J2 = ^jJ2{yo), J2 = ^^, Y2 = ^-Y2{yo), Y2 = ^^ (17) 
where J2 and Y2 are Bessel functions of order 2. The arguments yo and i/o are defined as: 



M / 4M2sin(x)2 , ^ 

^° = ^^^' ^° = ^V'^9cos(x)W3' (18) 



and coefficients M and K are defined as: 



, . 1^ -r. - 9cos(x)^°/3 

k y 4M2sin(x)2 ^ ^ 

Note that the expression for ( IT5|) is a function of x and M, only! This allows us to 
fully explore the analytic behaviour of the probability density which we shall do in the 
following subsections. In particular in the thin wall limit {x — > 0) A"^ takes the form: 

^"^ " n^M[Jf{M) + Y^{M)] ■ 
For the sake of completeness we also present the expression for the ratio: 

(21) 



and bm = 1/ ^/l + {a„JbmY. 

Note that for the thick wall x > 1 the expression for the probability density of the wave 
function varies significantly over the range of the wall. In this case the probability density, 
averaged over the domain of the wall thickness, may be more appropriate: 

m' = AlFm^ Fm = \{l + "^^). (22) 

Note that is a mild function of x, and while ^ 1 as x ^ 0, for x — ^ it approaches 
the asymptotic value \. 



A. Small and large Kaluza Klein Energy Expansion 

The expression ( IT5l) can be readily expanded in different regimes of values of x and M 
parameters. 

The expression (fT5l) has a universal behaviour: 

^ - , as M ^ 00 . (23) 

vr 

On the other hand one can also explicitly expand (ITS!) in the limit of small M. We have 
done so up to the fourth order in M expansion: 

Al^ = Mao{l + M'^ai) + 0{M^), for M<1 (24) 

where coefficients and ai are the following explicit functions of x: 

9cos(x)^°/3 

2(cos(x)'' + 2xsin(x) + 2cos(x))"' P2 



where 



2 cosixY^^ 

Pi = cos(x) ( 8 - 34cos(x)^ - cos(x)''(l + 541og( —r^)) 

\ Mexp(7) 

+2 X sin(x) (4-9 cos(x)^ - 4 cos(x)^) , (26) 

P2 = 18 cos(x)^°/^(cos(x)^ + 2xsin(x) + 2 cos(x)) . (27) 
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FIG. 1: Probability density A"^ plotted as a function of M for a small value of x = 0.1. 



The thin wall limit {x 0) produces the following values of ao and ai: 



"0 = - , "1 



i + log 



• AI exp(7) 



) , as X — > . 



(28) 



Note ( 1271) is a valid expansion for x < |, however in the thick wall limit (x |), Am has 
further suppressions for small values of M : 



"0 



9 (7T _ 3,U0/3 



27r2V2 



^1 ^ 9 V 2 



10/3 



as X 



(29) 



and thus Al ^M^. 

In figure [1] and [2] is plotted for a small and an intermediate value of x, respectively. 
For the sake of completeness we also give the expansion: 



^ = m2/?o ( 1 + /?i ) + OiM^) , for M < 1 

bm 



(30) 



where: 



where 



SttcosCx')^/^ 



4(cos(x)'^ + 2x sin(x) + 2 cos(x)) 



^^2 



(31) 



ri = 4x sin(x) cos(x)(2 — 11 cos(x)^) + cos(x)^(28 — 92 cos(x)^ — 17cos(x)'^) 

,2cos('x)2/^ 



4x^ sin(x)^(5 - 8 cos(x)^) - 108 cos(x)^ log( 



Mexp(7) 



r2 = 72 cos(x)"'^'^/'^(cos(x)'^ + 2x sin(x) + 2 cos(x)) 



(32) 
(33) 
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FIG. 2: Probability density plotted function of M for x = 1. 



whose thin and thick limits have the following form: 



4 ' 



/3i = -| + ilo, 



' M exp(7) N 



as X — >■ 0; 



4V 2 



X] 



Pi 



72 V 2 ' 



as X 



2 • 



(34) 



B. Resonances 



Another interesting behaviour of Am takes place in a particular window for a combination 
of parameter a; > 1 and small values of M: 



yo 



M I AM^ sin(x)2 



cos(a;)^/^ 



9cos(a;)io/3 



(35) 



or 



M Mtanfa;) 

< 1, > 1 



cos(a;)^/^ ' cos(a;)^/^ ^ ^ 

The probability density (|T5ll exhibits a resonance behaviour 3]. The resonance then appears 
when ^25' — i^'C* ~ 0, which for the regime (1361) takes the form: 



S 3 



(37) 



For the above expression we have employed the small argument ?/o ^ 1 expansion of Y2 = 
^ ^ Oivl'"^) and = ^ + C(?/o~^^^)- Eq. ([37D is satisfied approximately when 



3/2 
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yo~i(2n+l) (n = 1,2,---) or: 



M„~3cos(2;)^/3^n(n + l), n = 1, 2, ■ ■ • . (38) 

Note that the above relation for the position of resonances becomes less valid as n increases. 
When the condition (137|) is satisfied, the probability density (fTS!) is highly peaked at 

,2 18 18cos(x)i°/3 , , 

Al-^,= JJ, »1- (39) 



7r2 yl TT- 



To obtain the above expression, we set in f|T5|) : l^'S' — l^'C ~ 0, and employed the small Uq 

argument expansion of of J2 = ^^i^ ^ + C>{yl^'^) and = ^^2° 1" ^(l/o^^) which yields 

_ _ _ _ 5/2 

(along with (JSTD): Js^ - ~ — 3\72 ' "^^^^ the sharp fall-off of the density amplitude with 
the increase of the discretely valued M„ (138|1 . 

One can also estimate the approximate width AM„ of the resonances, by determining 

5/2 

the change in the parameter M (by |AM„) when Y2S — Y2C ~ J25' — J2C ~ using 
locally linear approximation. This condition yields: 



Mn 9(yo-|^) 3C0S(X)5/3V 4n(n+l) " 



Note that in this case the estimate for ^Mn, which effectively parameterizes the integral 
of the sharply peaked density probability over its width, takes the form: 

and is thus close to being M„ independent! Of course, this is a valid approximation only for 
the first few resonances whose width is still small, as we used only locally linear approxima- 
tion for the quantity Y2S — Y2C near M„. 

The most distinct and pronounced resonances appear in the range x ~ {1.35, 1.55}. In 
figure [3] the density probability as a function of M is plotted for the regimes of x where 
the resonances are pronounced. We have checked explicitly that our analytic results for the 
location, the amplitude and the width of resonances are in good agreement with numeric 
ones, as it is evident from table [B 



IV. MODIFICATION OF NEWTON'S LAW 

The analytic form of the probability density in turn allows us study explicitly deviation 
of Newton's law in the presence of the continuum states. The (bound state) graviton wave 
function (fTOj) determines the four- dimensional Newton's law, while continuum Kaluza Klein 
fluctuation modes (fT3|) are responsible for corrections to Newton's law. The four- dimensional 
gravitational potential between two point-like particles with Mi and M2, a distance r apart, 
can be written in the form (see e.g., jol): 

V = G„^(l + .), (42) 
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FIG. 3: Probability density plotted as a function of M for the value of x = 1.5. It exhibits 
distinct resonances. 





n 


1 


2 


3 


numeric 




0.05141 


0.08889 


0.12543 


analytic 


Mn 


0.05133 


0.08891 


0.12574 


numeric 


42 


680 


50.07 


10.12 


analytic 


42 

- '-Ill 


719 


48.00 


8.50 


numeric 




1.0710-5 


1.4710-^ 


7.1610-^ 


analytic 




1.0410-5 


1.5610-^ 


8.75 10-"^ 



TABLE I: The comparison of the numeric and analytic results for the location M, the amplitude 
and the width AM of the resonances for x = 1.5. 

where Gn — ^pianck^^Q correction to Newton's law is of the form: 

r e-^'W^^dm. (43) 
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Note that since NQ/k (eq. (fTT!) ) is only a function of x and \'ipm\^ (eg. (122!) ) is only a function 
of X and M = m/fc, the correction to Newton's law can be rewritten as an integral over M: 

'5=4/ e(-^'='-)|V^„|2dM, (44) 

and it is thus only a function of kr and x\ 

Note again that the averaged wave function density iV'mP (eq. fl22|) ) is of the form A'^Fm 
where Fm is a mild function of x and M, with — >■ 1 and Fm — '•l/2asa;— >0 and a; — *• |, 
respectively. In the following we shall focus only on factor and its effect on the structure 
of 6. 

The amplitude (eq. ffTSl) ) has a small M expansion of the form and fl27|) and 
it asymptotes to tt~^ as M — > cxd ((eq J2^ ). For a; < 1 it is a reasonable approximation 
to split the integral flH|) into two intervals M = {0, Mq} and M = {Mc, oo} where the 
form of the small M and large M expansion is used, respectively. A reasonable choice 
Mc ~ (vrao)"^, produces an approximate result which is in good agreement with numerical 
one: 

' = W^ijBy^'-^'-'''''''^(^okr + l)] 

+ -^[6-e(-*'^'='"V(Mcfcr)3 + 3(Mc;A;r)2 + 6Mcfcr + 6)] + ^e("*^^'=^H. (45) 
[krp nkr 

Again, note for r ^ {kMc)~^ the leading correction is of the form: 

k ap Goi 
^=NlJkrY^^^JkrY^- 

The intriguing possibility takes place for x > 1 where the resonance behaviour takes 
place. Note that in this case the resonances are highly peaked, and the integral over M 
can be approximated by a sum of expressions kN^"^ y4^AM„e~^"'^^ over n. Employing the 
analytic expressions from subsection IIII Bl we obtain the following correction to Newton's 
law: 



2 



^~zy:A/i + i;i7^e--"-, (47) 



'Mn kr 



TT V 4n(n.+l) 

n 

where M„ are locations of resonances (l38il and the approximately constant prefactor is 
obtained by using the limit a; ^ f in ([11]): kN^'^ -> ^cos^xf/^ and dH]): (A^AM„) ~ 

cos(a:)^/^-^l + ^n(ji+i) - Note that the prefactor is approximately constant and of order 
one. 

-2 



M„'s fl38|) are typically in the range of 10 and thus for a choice of also small k, and the 
range of distanced r ~ {Mnk)-^ the leading resonances can potentially contribute a sizable 
effect! In particular when choosing k ~ 10^ GeV, at distances that r ~ lO-^GeV"^ the thick 
wall limit with resonances M„ ~ 10~^ would produce an order one modification of Newton's 
law, while for the walls with x <1 the corrections ( H5f^6l) would be negligible, i.e. (9(10~^). 



V. CONCLUSIONS 

We have presented the first explicit model of a finite thickness domain wall, interpolat- 
ing between Z2 symmetric five-dimensional anti-deSitter vacua, where the graviton wave 
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function fluctuations can be studied explicitly for any thickness of the wall, parameterized 
by a parameter x = {0, |}. This allows us to determine explicitly the probability density 
of the Kaluza Klein fluctuations both for the small and large values of the Kaluza Klein 
energy. Notably, for x > 1 resonance behaviour emerges, which is most pronounced in the 
range x = {1.35, 1.55} and can be analysed explicitly. For a specific range of Kaluza Klein 
momenta (and anti-deSitter cosmological constant A) this resonance behaviour can signifi- 
cantly modify four- dimensional Newton's law by an effect of order 10^ larger than that of 
thin walls. 

While the concrete model is based on a specific, analytically solvable Schrodinger potential 
for the graviton fluctuation modes, we have not addressed the origin of the supergravity 
Lagrangian that would result in the BPS domain wall solution whose metric leads to the 
proposed Schrodinger potential. Certainly, this is an important outstanding issue. 

The analytic solvability of the model for the graviton fluctuations lends itself to the further 
study of fermionic and spin-one fluctuation modes. In addition, while we chose to study the 
fluctuation modes for co-dimension one (domain wall) configuration in five- dimensions, the 
study to other dimensions is readily available and is relegated to future work. 
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